Brightening the (130 GeV) Gamma- Ray Line 



Zhaofeng Kang, 1 '* Tianjun Li, 1,2 '^ Jinmian Li, 1 '* and Yandong Liu 1 '^ 

1 State Key Laboratory of Theoretical Physics and Kavli Institute for 
Theoretical Physics China (KITPC), Institute of Theoretical Physics, 
Chinese Academy of Sciences, Beijing 100190, P. R. China 
2 George P. and Cynthia W. Mitchell Institute for Fundamental Physics and Astronomy, 
Texas A&M University, College Station, TX 77843, USA 
(Dated: March 8, 2013) 

Abstract 

The gamma-ray line from dark matter (DM) annihilation is too weak to observe, but its obser- 
vation will uncover much information, e.g., the DM mass and an anomalously large annihilation 
rate ~ 0.1 pb into di-photon. In this work, we construct a minimal effective theory (EFT) incor- 
porating DM and heavier charged particles. A large annihilation rate is obtained from operator 
coefficients with resonance or strong coupling enhancement. The EFT is stringently constrained by 
the XENON100 and WMAP data. Without resonance, Dirac DM or colored charged particles are 
ruled out. It is pointed out that the di-gluon mode may correctly determine the DM relic density. 
Interestingly, this framework also provides an origin for the Higgs di-photon excess at the LHC. 
We apply the general analysis to the NMSSM, which can elegantly interpret the tentative 130 GeV 
gamma-ray line. A top-window model is also proposed to explain the gamma-ray line. 
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I. INTRODUCTION AND MOTIVATIONS 



The existence of dark matter (DM) has been confirmed by its gravitational effects, and its 
energy fraction ~ 25% today is also measured. However, the conclusive evidences that may 
reveal the DM particle properties are still absent. Among a variety of (indirect) detecting 
objects on DM, the gamma-ray from the DM dense region (such as the center of the Galaxy) 
is especially promising by virtue of weak astrophysical influence on its propagation [1]. Of 
particular interest is the monochromatic gamma-ray line, which has rather clear background. 
But it is highly suppressed because the DM \ can only annihilate to photons via the charged 
loop. 

However, once such a spectral line is observed, it will uncover very important information 
of DM. In this article, we assume an extracted DM mass from E y and an anomalously large 
annihilation rate (crv)2j ~ 0.1 pb into di-photon (it is taken as a referred value throughout 
the work, unless specified), then attempt to reconstruct the DM properties and dynmiacs 
to the most extent. Inspired by the recent discovery of a gamma-ray line at E y ~ 130 
GeV, which is claimed in the Ref. [2, 3] after re- analyzing the Fermi Large Area Telescope 
(FERMI-LAT) data published in 2009 [4], it is conjectured that the line may originate from 
DM annihilating into gamma. Best fit of the data shows a DM of mass around 130 GeV and 
annihilation rate at level 0.1 pb. Later independent analysis also confirms the line [5]. The 
line has a sharper peak which is hard to explain by FERMI-bubbles [6], while DM+DM— > 77 
gives a better fit [5]. The Ref. [7] also shows a strong evidence of the gamma-ray from the 
inner galaxy and draws a similar conclusion. This line has received much attention from 
astrophysics [5, 6, 8, 9] and particle physics [10]. In spite of queries [8], the gamma-ray line 
from DM activity itself is of great theoretical interest, and deserving a deep study. 

The topic can be studied in the effective theory (EFT) framework, by minimally including 
an operator acy^xC^C where C is the charged particle. The anomalously bright gamma-ray 
line is due to large ac- We further demand the EFT be compatible with other constraints 
on the DM, i.e., the WMAP and XENON100 bound [12]. Independent of the mechanism 
generating large ac, we can arrive: 

• The charged particle C in the loop should be heavier than the DM, otherwise it would 
render too large annihilation rate into CC, which leads to too small DM relic density. 
On top of that, the injection from such a large flux of charged particles into the 
cosmic-ray probably has been excluded by the PAMELA. 

• The charged particle carrying both QED and QCD charges needs careful inspections. 
Along with the di-photon annihilating mode, there is an enhanced di-gluon mode 
with estimated rate (<jv)2g ~ 0.1(aj./a 2 )(av)2-y — 1 pb, which makes an illustrative 
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coincidence. 



The large ac can be generated through Breit-Weigner resonance mechanism, or the 
strong interaction between DM and the charged loop. For the former scenario, proper- 
ties of the scalar/vector resonance can be further stringently restricted by symmetries, e.g. 
the CP, and the above consideration. The next-to-minimal supersymmetric standard model 
(NMSSM) [13] is a good realization of this scenario. We find it is capable of interpreting the 
tentative 130 GeV gamma-ray line. For the latter scenario, the XENON100 bound excludes 
the Dirac DM, as well as both Dirac and Majorana DM if the charged paticle C carries 
color. Interestingly, in any scenario, the possible SM-like Higgs h to di-photon excess at the 
LHC [14] may share the same origin, if we incorporate the operator a^hC^C. 

This paper is organized as following: In the section II, we perform a general analysis based 
on the minimal EFT. In the next two sections exploration on the enhancement mechanism 
is presented. The Section V includes the conclusion and discussions. And some necessary 
complementarity is casted in the Appendix. 

II. GENERALITY AND GUIDANCE 

As is well known, the DM can not directly annihilate into photons due to its QED 
charge neutrality, while transition at the loop level is generically highly suppressed. So, it 
is nontrivial to obtain an abnormally large annihilating rate, saying {av)^ ~ 0.1 pb. Some 
more powerful model independent statements can be made, if it is further combined with 
other aspects of DM. To see that, we consider the minimal effective operators [42] relevant 
to the gamma-ray line anomaly 

FermionDM: acxTx^TC, gqxTx&C, a wX TxW + W~, (1) 
Scalar DM : a t ^xCVC, aw^xW+W, (2) 

where the gamma matrix Y C {1, 7 5 , 7 M , a^ u }. Lorentz and SU(3)c x ^(1)qed x CP sym- 
metries are implied, and some operators will vanish due to these symmetries. C is a charged 
fermion and C is a charged scalar, both of which are not confined to the SM. However, in 
the sense of inducing DM annihilating into gamma, the scalar loop is not as effective as the 
fermionic loop, unless there is a large enhancement from highly charged particles such as 
a double charged scalar. Note that operators containing H + W~ where H + is the charged 
Higgs from 2HDM-like model are not included, since their contribution are always putative 
null, e.g., in the non-linear unitary gauge the vertex H + W~ , ~f vanishes [15]. 

To achieve a large (<jv) 2i and maintain the main merit of DM dynamics, the candidates 
running in the charged loop are more or less selected. Denoting by Si the set of DM 2—7-2 
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annihilation mode and without loss of generality, let Si be the one from which di-photons 
come after closing the charged states to form a loop. Some cases arise: 

• If Si is the on-shell type with final states XiXi, then (uv)s 1 ~ 10 4: (av)s 1 , 2l ~ 10 3 
pb. Injection from such a large flux charged particles would have been observed by 
PAMELA [16] from the significant excess of the positron or anti-proton flux. On top 
of that, the DM relic density would be too small, unless we consider the subtle thermal 
Breit-Weigner enhancement effect which is active only today [17]. 

• If Si is properly off-shell, then the above problem is resolved since the Si is forbidden 
today. But we have to examine its annihilation rate at the early universe, i.e., compar- 
ing the rate of the forbidden annihilation mode {(tv)s 1 ,t j with 1 pb, where Tf = m x /xf 
with Xf ~ 25 the typical decoupling temperature of DM. Generically one can expand 
the annihilation rate as [18] 
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i = (a + b/xf)v 2 , (3) 



where the final two-body phase space gives the velocity of out-going particles in the 
CM frame, v 2 = (1 - z 2 + z 2 v 2 el /4) 1/2 with z = mx 1 /m x and v re \ the relative velocity 
of initial particles. For a properly large z > 1, the relic density of DM is [18] 

__, 1.07 xlO 9 *, z a 

nh = TB > J ~ a e y > 4 

gl /2 M Pl J fi-Xf 

where /x_ = (1 — 1/z 2 ) 1 ^ 2 . The relic density is very sensitive to z, e.g., from z = 1.05 
to z — 1.10, it increases roughly one order. But in principle it is possible to obtain 
{<jv) 2i ~ 0.1 pb and {cv)x 1 x 1 — 1 pb simultaneously, if we accept large fine-tuning. 

If z is large enough then the forbidden channel is completely ignorable. As a con- 
sequence, we need a new (dominant) channel S 2 to reduce the DM number density. 
In actual model building, it naturally happens. But a more interesting case arises as 
following. 

If the charged particle also carry color charge, then the 27 mode is subdominant to 
the two gluon mode: 

MS ~ 0.l4 = 0(20), (5) 

which is estimated in light of fermions with unit charge, and the origin of numerical 
factor 0.1 can be traced back to the property of charged particles. This numerical 
coincidence means that if the two-gamma rate is ~ 0.1 pb, the right relic density is 
achieved via the di-gluon mode. For the proof of that the di-Z mode is at most the 
same order of di-photon mode and thus irrelevant, see the Appendix A for details. 
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Generically, we will have several annihilation modes producing gamma line, \X ~~ >* Hi lX 
with X = Z, h. The first mode creates a line at E 1 = m x while the second mode produces 
a line with lower energy E' « Ey(l — mj^/m 2 ,) 1 ^ 2 . They have comparable cross sections 
except for a very significant phase space suppress. Therefore we focus on the two gamma 
final states, since the two modes quantitatively differ only by some constant, such as the 
difference between e and #2- However, how to distinguish the two lines is very interesting as 
discussed in Ref. [19]. 

To end up the general discussion in EFT, we would like to mention that there is a possible 
relation between the di-phonon excess for the Higgs search at the CMS/ATLAS [14] and for 
the DM search at the sky. The common point is the new charged loop. Through the same 
loop alone which the DM annihilates into two photons, the SM-like Higgs h can decay into 
two photons with appreciable width, if the coupling to the Higgs is significantly. It is can be 
described simply by further including the effective operators a^hCC or a^hC^C. However 
we are not going to discuss this in detail due to its triviality in the EFT. 



III. ENHANCEMENT FROM RESONANCE 



In this section, we present the effective analysis by specifying the role of resonance. The 
simple top window model is constructed, and in particular we survey its implication on 
conventional supersymmetric model such as the NMSSM. 



A. Scalar resonance 



We consider the scalar resonance which appears almost everywhere in models with ex- 
tended Higgs sector. To get a sufficiently large enhancement, s— channel resonant anni- 
hilation is the most conventional mechanism. In this case, the resonance <p takes mass 
m^ ~ 2 777. rjM> an d the cross section manifest of the enhancement can be parameterized as 

T i T F 1 , , , l2 « 2m Y 1 

where r = 4m 2 , jm\ and 7 = (T^/m^) 2 <C I. Here a stands for an effective coupling and 
will be specified in concrete examples, while Tj t p takes 1/2 or 1/4 and so on, standing for 
the average of initial degree of freedoms or the symmetry factor of final states. 

It is convenient to define fg = 1/((1 — r) 2 + 7), then the DM annihilating cross section 
into XiXi (thought 0) can be rewritten as 

(av) Si = I^^L\M(XX -> 0)| 2 1X0 -> XiXi). (7) 
ml m x 
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Without loss of generality, we take X 2 as the dominant mode, and immediately get the 
upper bound of the branching ratio of <fi decay (to particles other than di-photon): 

Br(0->X 2 X 2 ) _ (av) S2 
Br(0 -> 77) (H27 ~ ' 

To arrive it we have set (<jv)s 2 ~ 1 pb as the standard annihilation rate as well as the 
referred value (crw) 27 ~ 0.1 pb. Therefore we get a model independent bound Br(0 — > 77) > 
10%, as negates the resonance from simple two-Higgs-doublet-model (2HDM) by virtue 
of their considerably coupling to fermions or light massive vector boson (it is absence for 
CP-odd Higgs). This bound has far-reaching implication on the collider. Provided that 
the production cross section of <fi is sufficiently large (for example, when the charged loop 
meanwhile carries color as in the top- window model discussed later), the gamma-ray line 
observed at the sky predicts a clear di-photon excess at the peak around 2m Y (^ 260 GeV) 
at the LHC. 



1. Effective analysis 

In light of previous arguments, we need some rather heavy charged particles, while in the 
SM top quark and W boson are the only two charged particles of mass around the weak 
scale. Accordingly, a m x > m t hints a new charged particle. When m t > m x > mw, the 
top quark will open a unique window. Otherwise, W may run in the charged loop provided 
a vertex <j)W + W~ , that implies nontrivially the identity of 0. 

Since the effective operators listed in the previous section are ascribed to the integrating 
out s— channel scalar resonance, the set of possible operators can be reduced greatly: 

Fermionic DM: a cX (^)xC(l 5 )C, a d x(j 5 )x&C, a w xxW + W~, (9) 
Scalar DM : a c ^xC(j 5 )C, a dX j xC j C, a w x ] xW + W~. (10) 

The coefficients a's are proportional to f l J 2 /m 2 ^. And 7 s in the parenthesis may or may 
not appear, depending on the CP quantum number of <f. But specified to fermionic DM, 
7 5 must be inserted so as to make the present DM annihilating rate avoid acute velocity 
suppressing, v 2 ~ 10 -6 . This means the must be CP-odd, denoted as <$a hereafter. As an 
immediate consequence, the W— window is closed. Additionally, C'C should be understood 
as C^Cr, under CP transformation C L , R —> Cr/l- However, QED does not change the 
chirality, thus we need a further large LR mixing. As an example, the stop-system in the 
SUSY just satisfies those requirements. On the contrary, for scalar DM, no matter complex 
or real, the <fi should be CP-even and denoted as <ph- Hence the 7 s should be removed. 
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We would like to add some further remarks. Firstly, the Lorentz and SM-gauge invariance 
force either transforms non-trivially under the SU(2)l symmetry or mixes with such states. 
Secondly, the constraint on the interactions indicated by the Eq. (8) should be satisfied. 
Finally, the also mediates the tree-level DM-nucleon interaction in the presence of a 
top window. Although this contribution is suppressed by velocity for the fermionic DM, 
the scalar DM requires inspection [20]. The resulting DM-proton inelastic scattering cross 
section is 

r2 n 2 p a t m p 
asi = — f P> ^27^2^? 

where fi p ~ m p is the DM-proton reduced mass and ~ 0.83 [20]. Note the enhancement 

1 12 

factor f B in a t is removed when we are calculating the DM-nucleon recoil rate using the 
effective operator. The present exclusion on f p is 10 -8 GeV -2 for DM of mass 100 GeV, put 
by the XENON100 [12]. It implies the upper bound 

a, < 3 .2f pmx m t/mp = 2.8 x l«r» (jjp^) (lo^v) &V "' < 12 > 

which places a rather strong constraint for the top-window model. 

Having outlined the most essential profile of WIMP that potentially has bright gamma- 
ray lines, we continue to make some quantitative discussion. Effectively, through the charged 
loop, Lorentz and CP invariance leads to the following operators for the CP-even and CP-odd 
respectively 

a^^F^, a ^±^r (13) 
4vr 4Ai M 4tt 8A 2 m V ; 

with a ~ 1/137. h^cc is the coupling constant between and charged particles. Factoring 
out the loop factor and couplings, moreover multiplying 1/4(8) for later convenience, the 
A12 can be much below the weak scale (it is even enhanced by color or electric charge). 
Concrete expressions for the effective scale are casted in the Appendix A. 

Now we are at the position to evaluate the DM annihilation rate into gamma pair. De- 
noting by the Feynman rules (see Fig. 5 for label) of Eq. (13), allowing for off-shell 0, 
they are respectively given by [43] 

CP-even: r^( Pl ,p 2 , P) = \pi • T*gT ~ (j>1vl , (14) 

CP - odd : I^fa, q 2 , P) = ^^^e^{p x ) a {p 2 ) p . (15) 

After averaging initial states and summing over out-going photon polarization states, the 
resulting real scalar DM annihilation rate is 

] ^-l^l^^^^ x ~ 0m {500){-Tr) { 0.3 ) Pb ' (16) 
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where gu xx is the coupling constant between <fi and DM. Owing to the scalar DM, we have 
parameterized the massive coupling as g^ xx —> 2g' l j > m x . For the Majorana DM, 




FIG. 1: Dark matter annihilates into 27 via top loop. The cut denotes for annihilates into top 
pair. 

2. The top-window model 

In the standard model (SM), lying on the top of the fermion mass ascending order, the top 
quark may hide some new dynamics. It is thus interesting to conjecture that the dark sector 
may have a close relation with the top quark, i.e., the dark sector only (strongly) interacting 
with top quark in the fermion sector [44]. Then it is reasonable to expect enhanced gamma- 
ray line via top-loop at levels close to the present experimental sensitivity. And interesting, 
the latest tentative gamma-ray line from dark matter annihilation [3] requires a DM of mass 
about 130 GeV, just lies within the top-window. Therefore, as a concrete example of the 
general effective operator analysis, we focus on a top-window model equipped with scalar 
resonance. 

Due to the chiral structure of the SM, it is most likely that the <fr dwells in an extra Higgs 
doublet H' with hypercharge +1/2 so as to couple the top quarks at dimension-four level. 
A simple effective model, which might be a remnant of top dynamics reads 

C = l - X 2 {(% + XH'H + c.c. + Kl \Hf + k 2 \H\ 2 ) + V(H, H') + y 3 Q 3 H't R . (18) 

where H is the ordinary Higgs doublet developing a vacuum expected value (VEV) v = 
174 GeV. By contrast, H' should have a VEV (H') <C v which makes the term (pW + W~ 
neglectable [45] . Here we assume that the VEV hierarchy is realized by well organized Higgs 
potential V(H,H') as discussed in [20]. In addition, since we have ascribed the unique 
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significant coupling H'Q^tn to some unknown dynamics, we are free of the FCNC-problem 
induced by family changing Yukawa couplings H' 'JlIr- in summary, after the EW-breaking, 
we are left with the relevant terms 

C D i (n\ + k 2 v 2 ) x 2 + ^j={\v) X 2 <P + -^={n 2 v) X 2 h + ^ X 2 h 2 + \^J=^ L t R + c.c)j . (19) 

In light of the general analysis, the scalar DM the <j> is identified with the CP-even component 
of the neutral boson in H 1 . If we drop terms involving H', the model just recovers the Higgs- 
porting model (see a recent discussion [21]). But since k 2 is irrelevant on our purpose, so 
we can turn it off to reduce parameters. 

Now we turn our attention to the phenomenological aspects of the model. First of all, 
the di-photon rate, in terms of the parameterization in the Eq. (17), is determined by 

, A V r- 

94>xx = ^3/2 —> h ^c = 2/3/ V2. (20) 

Taking A ~ y% ~ 1 and = 500 leads to a rate ~0.07 pb. And the XENON100 constraint 
Eq. (12) has been arranged to be satisfied by means of rather large fs thus smaller couplings. 
Next, the mass difference between m t and m x is at a few ten percents level, so the forbidden 
annihilation mode (to top quark) is ignorable in terms of the Eq. (4). The di-gluon mode 
may properly account for the relic density. From Eq. (A10) and Eq. (A7) it is estimated 
that 

(av) 2l a 2 4A c 2 Qt ' 1 ' 

with Q t = 2/3 the top quark charge. 

An alternative model is to assume the dark sector consists of a SM vector-like Dirac 
fermion pair (ip,ip c ), which is the the SU(2)l doublet and carries hypercharge ±1/2 respec- 
tively. Extra singlet fermion S is introduced, then the model is 

C = (XS^H 1 + m^ c + M S S 2 ) + V(H, H') + y s Q 3 H't R . (22) 

The fermonic DM brings important difference. The (f> should be the CP-odd component of 
H', which does not couple to W + W~ thus allows a large (H') (but we have to ensure small 
H' — H mixing, otherwise (fbb is too large). It leads to a phenomenologically viable singlet- 
doublet mixing Majorana DM. It can be regarded as an effective model of the supersymetric 
model studied as the following (but top loop will be replaced by chargino loop). 

B. The NMSSM 

In the supersymmetric standard models (SSM), in addition to the (Majorana) neutralino 
LSP dark matter of mass around m Z) a wealth of new heavy charged particles, CP-odd 
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Higgs and hence a viable resonance 0a, are furnished. It is thus of particular interest to 
investigate whether the above set up can be realized in the SUSY. Furthermore, as stressed 
at the beginning of this section, the coupling of the resonance to light states are stringently 
constrained. In spite of difficulties in the minimal-SSM, such an invisible <pA can be readily 
accommodated in the NMSSM following the singlet limit. While the examination on light 
dark matter limit has been done [22], our scenario has not been currently examined yet as 
far as we are aware. 
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FIG. 2: The LSP (with acceptable relic density 0.09-0.12) at the plane of (av )^ — Ogj- As one can 
see, the FERMI gamma-ray line can be accommodated even under stringent XENON100 exclusion. 
The colored notation is of the SM-like Higgs mass. We obtain the results by using the programme 
NMSSMtools [34]. 

Owing to the singlet sector, the NMSSM presents a clear realization of our scenario. For 
our purpose, the model is none other than the Z%— NMSSM: 

W D XSH u H d + ^S 3 , (23) 
-Csaft D m 2 s \S\ 2 + (\A x SH u H d + A K ^S 3 + c.c.) , (24) 

Consider a slice of the parameter space: (A) A ~ 1 moreover k < A and tan ~ 3 (favored 
by naturalness to enhance the SM-like Higgs mass [23]), further the v s = (S) gives rise to 
the Higgsino-like charginos (with mass roughly \i = Xv s ~ 200 GeV) which replaces the 
top quark in the gamma loop; (B) The 130 GeV LSP dark matter has significant (even if 
not dominant) singlet component, and kS 3 /3 provides the vertex nA s S 2 with unsuppressed 
coupling. (C) The highly singlet-like CP-odd Higgs A\ ~ A s of mass around 260 GeV offers 
a proper resonant enhancement. Consequently Abb and AhZ can be sufficiently suppressed, 
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and its only significant coupling is to the chargino ~ \A\H^H~^ . This hiding A\ is a key 
to reconcile the WMAP and FERMI, and the situation of this parameter space to interpret 
FERMI is shown in the Fig. 2. We have restrict falls into 255-265 GeV. As one can see, 
only a very small portion of the points pass all constraints, labeled as the FERMI-region. 

We close this section by making some further comments. The above parameter space is 
a portion of the natural NMSSM [23]. However, as observed there, the Higgsino generically 
occupies a large proportion of LSP, and its relic density is too small while (Xgj is too large. 
To circumvent those problems, we may have to tune the parameters to obtain a viable LSP, 
as is reflected in the Fig. 3. From model building, we may simply go to the singlet-port 
dark sector by simply adding r]S§ 2 /2, where the lighter $ component is the dark matter 
candidate. 



Distribution of (av) 2y 



Distribution of (ctv) 2y 




480 490 500 510 520 530 540 
M A (GeV) 



550 




270 



FIG. 3: The distribution of the di-photon rate on input coordinates, left: Ma — kA k ; Right: 
[i — Mi. The solid purple circles stand for parameter configurations satisfying WMAP, XENON100 
and FRRMI. They are well-tuned, since they are rather discrete even in the preferred window. 
Other parameter settings: A = 0.65, k : 0.13 — 0.16, tan/3 : 1.3 — 1.7, f^Q 3 = m t R = 1000 GeV, 
A t = while A R varies between -400 to -300 GeV. 



C. Vector resonance 

In this subsection we turn our attention to the vector resonance. There is only one 
neutral massive vector boson in the SM as well as its simple extension. Such a resonance 
means the DM mass should be around m^/2. However, in models with extended U(l)x 
local symmetries, new resonances with free masses are expected. Such models have been 
proposed in [24, 25], where spectral lines from 7Z and/or 7/1 final states are predicted. It is 
believed that the Landau- Yang theorem [26] excludes the di-photon mode. Nevertheless, this 
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theorem does not exclude the vertex Z77 with off-shell mz, please see the Ref. [27] discussing 
anomalous three vector boson coupling within SM. And a concrete evidence for such coupling 
can be found in the earlier calculation of neutralino annihilation to di-photon [15], where a 
pole from mz indeed exists. 

This pole has deep relation with the anomaly of the theory, and it is non-vanishing if 
and only if the axial-coupling between Z' and the charged massive fermions is present (mass 
splitting is needed to spoil thorough anomaly cancellation). As a consequence, in Eq. (9) 
we are left with the only one operator built from fermionic DM and charged particles: 

acXl^xCl'YC, (25) 

which indicates Z' can not come from vector-like theories. While for the complex scalar 
DM, the relevant effective vertex is x d^X^a-i which renders the DM annihilation suffering 
from velocity suppress. So we do not need to consider it here. 

Whatever the charged particles are, the effective Lagrangian can be built by Lorentz 
and QED invariance. This leads to the following dimension-six operators 

a d'zCC 7 i t-i rpua 771 , a 9ZCC 7/ H j?ua 771 




47T 4A2 



a 9ZCC 7/ i« t?vol 771 1 a 9ZCC 7/ i« 77W 77. \ (cya\ 

4T1at^ f Fa " + ^^ Z ^ F (26) 

where the symmetric and antisymmetric 2-rank tensors are defined as {Zs/aI^u — d^Z' y ± 
d u Z' fl . Under C and P symmetries, the vector field transforms as 

CVC~ X = -V, PV{x,t)P~ 1 = (-l)"V(-x,t), (27) 

and Pd^P^ 1 = (— 1) M <9 M . As a result only A 12 term conserves the CP, so in this work we 
only keep them (which is consistent with the presence of 7 s in Eq. (26)). From the effective 
Lagrangian, the annihilation rate is calculated to be (for illustrative purpose we only show 
the Ax-related part) 

2 2 / \ 2 

/ \ 29z'xx9z'cc fs ( m x\ 1 
{av)2 1 ^ a 



16tt 2 64tt V Ai / Aj 



1 



,o.o3(^) 2 (A)(_ a _ ) ^l^X)' pb , (28) 

On the other hand, from direct calculation similar to the Ref. [15], we get 

;27 ^ NlQt (gz> xx 9z<cc) 2 ^S lAl/2{T)l2 - (29) 

X 
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This leads to the effective scale 

A 1 = m^/2Q c N 1 c /2 \A 1 ^(r)\ 1 / 2 , (30) 
where r = mc/m x and A 1 ' 2 {t) = r(arctan l/\/r — l) 2 > 1 for r > 1. 

IV. CHARGED-LOOP PORTING DARK MATTER 

In this section we consider the scenario in the absence of resonance. Viewing from the UV- 
completion level (illustratively showed in Fig. 4, where a Z 2 symmetry can be consistently 
assigned on DM and charged particle), the operators list in the Eq. (9) are generated in the 
t— channel. To get a large annihilation cross section without resonance enhancement, it is 
expected that there is rather strong coupling between DM and the charged particles. That 
maybe consistent with the composite dark matter. However, in the case of decaying DM 
scenario [28], a large coupling constant can be avoided. 




FIG. 4: Scalar /Fermionic dark matter annihilates into 27 via a charged loop. Z2 symmetry can 
be assigned to the dark matter and charged particles. 



A. Annihilating scenario 

In all cases, the heavy charged loop is the major port between DM and visible sector. 
This is a reminiscence of the dipole dark matter theory, where the DM-photon interactions 
are the leading order of DM-visible interactions. Considering the fermonic DM case, the 
effective operators up to dimension- seven should be incorporated: 

- y (*VX) F>iU , "4 (X<V7 5 X) F»\ (31) 
■^XXF^F^ iJL^xF^F^, (32) 
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Eq. (31) are dubbed magnetic momentum DM (MDM) and electronic momentum DM 
(EDM) respectively, with \ x /d x the magnetic/electric diploe momentum. If the DM is 
a Majorana fermion, Eq. (31) vanishes and we only have to consider Eq. (32). Otherwise, 
we expect all operators listed in Eq. (31) and Eq. (32) are comparably important, since they 
are generated at the same loop-level. The naive dimension counting gives 

1/A ^16^i' As-AW^/A,)^, (33) 

up to an overall operator coefficients loop factors /(m^/m^.), whose exact expression is very 
involved due to multi propagators, and we leave it for further work. 

Before dealing with the annihilating rate into gamma pair, we consider the possible bound 
on operators. In spite of loop suppression, the M(E)DM has long-distance interactions, 
mediated by the photon, which lead to a great enhancement on ogi- Especially, the EDM 
has a further 1/v 2 ~ 10 6 enhancement [35, 36], consequently the XENON100 tightly bounds 
on them. The DM-proton (only for proton by virtue of the QED mediator) cross section is 
°si = a ^x~ ^ n ^ e light °f R e f- [36], for a 100 GeV DM, the upper bound is roughly 

A x < 10~ 19 ecm « 1.5 x 10^ 6 GeV" 1 , (34) 

which indicates ^c/fi'Lc ~ 1 TeV. In a way similar to the one given in Section. IIIAl, the 
correlation between operator coefficients leads to an estimation rate of the XX ~* 77 process 



Mfrr-^AJ £ 10~ 6 pb, (35) 

which is far below the sensitive bound. Thus the DM can not be a Dirac particle. 

Direct detection possibly gives a second constraint, no matter Dirac or Majorana DM. 
If the charged particles also carry color charge, then after replacing the QED field-strength 
with the gluon field-strength, we get the dimension-seven operators 

fc^XxG^G% fafi^X^xG^G^, (36) 

Again from naive dimension estimation, /g j5 ~ ]q ~ 167r 2 /eA|. Unlike the MDM or EDM 
operators, they always lead to direct detection signals. Concretely, the first operator gives 
contribution to Ogj in the form of Eq. (11), with [29] 

— * -hafto < ID" 9 GeV- 2 /™n, (37) 

which implies the lower bound A3 < (IG^fc/e) 1 ^ 3 ~ 3.7 TeV renders again a very small 
(0^)27 ~ 10 -6 pb. Now we can draw the conclusion: in the charged loop ported DM scenario, 
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if we expect a large annihilation rate to gamma, XENON100 excludes both a Dirac DM and 
colored charged loops. 

We would like to comment on the scalar DM case, the (S^d^S) d u F» v leads to contact 
interaction between DM and nucleon and thus safe. In fact, such model [11] has been recently 
proposed to explain the anomaly using the scalar DM. 



B. Decaying scenario 



The scenarios discussed previously involve either some tuning or strong couplings, while 
decaying DM gives an alternative more natural choice. We close the paper by giving a short 
comment on this scenario to explain the FERMI gamma-ray line. Here, the scalar DM is 
more or less favored. At the two-body decay level (three-body or more leads to too wide 
spectrum to account for the peak), the scalar DM can decay into the gamma pair while the 
fermionic DM can not. The unique effective operator is written as 

~xhF^, (38) 

but this time the \ obtains a TeV scale VEV v x which breaks the Z 2 symmetry and leads to 
the scalar DM (the real part of x) decaying to gamma pair. To fit the data, besides a mass 
of DM should be around 260 GeV, we further need its extremely narrow branching decay 
width to gamma pair ~ 10 _29 s _1 [28]. The small decay width can be achieved by lifting 
the mass of charged particles running in the loop, typically A c — > Mqut (see an example 
in [33]): 




3 x 10 14 GeV 



y*- 1 . (39) 



A J 

Some comments are in orders. Firstly, in our notation the A lies below Mqut about two 
orders, but actually it may be compensated by loop factors and large Yuawa couplings (set 
to unit in the above estimation). Secondly, the relic density is a generic problem for decaying 
DM, but non-thermal production such as Ref. [30] using freeze-in mechanism [31] may offer 
a solution. Last but never the least, the 130 GeV gamma-ray line in decaying DM scenario 
is easily compatible with the sharp excess in PAMELA positron fraction [16], which can be 
interpreted by a leptonic decaying DM (saying to e + e~) with lifetime r ~ 10 26 s and mass 
around 200 GeV. In model building, one has to introduce relevant dimension-six operators 
for DM decay to leptons and adjust the branching ratio to fit both datas. 
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For the fermionic DM, the Lorentz invariance force the presence of second fermion (in 
the final state) in the operator of decay. As a case in point, in the (lepton number violating) 
R— parity violating SUSY, the gravitino has the following two body decay modes [32]: 

G^z/ + 7, W + +£~, Z° + u, (40) 

where gravitino mass is mg = 250 GeV and the branching ratios are respectively given 
by 0.03, 0.69 and 0.28. Although it fails to explain both PAMELA and FERMI gamma- 
ray line, the first mode can explain the latter given proper decay width. Of particular 
interesting, the neutrino (from the third mode) and gamma signal appear simultaneously, 
and the accompanied neutrino signal may be detected and thus provides a complementary 
detect method for this scenario. 

V. CONCLUSION AND DISCUSSION 

The gamma-ray line from DM annihilating in the galaxy center generically is well below 
the present detectable level. However, once the observation, from it we are able to extract 
very important information of the DM properties/dynamics. In this work, we present a 
minimal effective theory framework to understand the anomalously bright gamma-ray line 
from dark matter activity: 

• In the EFT large annihilation rate is ascribed to operator coefficients with resonant 
or strong coupling enhancement. 

• Due to the XENON100 bound, Dirac DM or colorful charged particles are ruled out 
in models with only strong couplings. 

• If the charged particle in the loop carry SU(3)c charge, the di-gluon annihilation mode 
is about one order larger than the di-photon mode, that may properly account for the 
relic density. 

• The SM-like Higgs may share the same charged loop, and therefore provide a source 
of Higgs di-photon excess at the LHC. 

Applying the general analysis to the NMSSM, that is proved to accommodate neutralino LSP 
with large annihilation rate into di-photon and interpret the tentative 130 GeV gamma-ray 
line. Top-window model is also proposed to explain it. 

Although not the central points of this work, we would like to end up by commenting its 
very promising collider detection prospect, if the 130 GeV gamma-ray line from DM activity 
will be confirmed. In light our general analysis in the text, at the LHC or Tavertron (but 
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beyond LEP) one can expect new light color-singlet charged particle C can be produced: 
qq — > C'C. While beyond the 130 GeV line and consider more wide scope, the LHC could 
put very strong exclusion on the model with colored loop. 
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Note added 

In the completion of this work, we note the appearance of work [40], they also note the 
importance of the di-gluon mode in the determining of DM relic density and more relevant 
phenomenologies are discussed there. The Ref. [41] specifically studies the 130 GeV gamma- 
ray line from LSP annihilation in the NMSSM, taking a quite similar idea to ours. We 
greatly thank the authors for sending us the updated version of NMSSMtools in which the 
CP-odd mass correctly adopts the running mass. Using it we re-scan our region, and find 
different results than [41]. 

Appendix A: Effective vertex from the charged loop 
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FIG. 5: fermion loop induced neutral scalar decay 
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1. -> 27 

In this Appendix we present the procedure of calculating the coefficients of the effective 
vertex used in the Section III A, quoted for convenience: 

4vr 4A/ h ^ ' 4tt 8A 2 ^ ^ ' { ' 

The calculation is similar to the case of the Higgs with general couplings. For definiteness, 
here we focus on the process depicted in the Fig. 5, where CP-even resonance <ph decays to 
gamma pair with single charged fermion running in the loop. At the one hand, the direct 
calculation of decay width gives 

1 & -+ TY) = 1 2N c Q 2 c h 4>cc A h 1/2 (r) 1 2 , ( A2) 

2507T J ' 

with N c = 3 the color factor and Qc the electronic charge of C. The loop function A(t) is 

A h 1/2 (r)= ^ [r + (r - l)/(r)] , (A3) 

Jarcsin 2 Jt t < 1 



i 

only depending of the ratio r = /4m 2 7 . In general, the <fth can be off shell and therefore 
the replacement — > P 2 , with P the four-momentum of <ph- For the CP-odd resonance 
tpAi the decay width takes the same form as the Eq. (A2) but the loop function is different: 

Af /2 = 2r-^f(r). (A5) 

On the other hand, from the effective operators in Eq. (Al) one can calculate 

a 2 h 2 a 2 h 2 
r <*^> = iiHi£sf<' T ^rt = w£k*'- (A6) 

So, comparing Eq. (A2) and Eq. (A6), we eventually get the effective suppressing scales 

A 1 = m x /\2N c Q 2 c A h 1/2 \, A 2 = m x /\2N c Q 2 c Af /2 \, (A7) 

where = 2m x has been used. Those expressions justify the statement that the Ai j2 can 
be much below the weak scale (or m x ). 

2. <p -> 2G, 2Z 

There are two other effective vertex require attention. First, if the charged particle also 
carries color charge, then the resonance can decay into two gluons. Repeating the procedure 



log^ 



VT- 
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FIG. 6: The loop function \A\ 2 varies as variable r = rn^/Arn^. Dashed line: for the CP-odd (j>\ 
Solid line: for CP-even cj). The former is always several times of the latter. 



dealing with the vertex (fihFF and <PaFF, the corresponding coefficients parameterized ex- 
actly as the Eq. (Al) except for the replacement a — > a s and A» — > Aj )S , we are able to get 
the precise effective operators. It is straightforward to get the gluonic partial decay width 
from direct and effective calculations: 



r 
n 



h ->■ gg) 
*>a -> gg) 



a l m 4> h 
32tt 3 

32vr 3 



^CC*^i/ 2 l 



128tt 3 A? " l hi 
128tt 3 A 3 " l A- 



Then we obtain the effective scales 



A liS = m x |/K 



1/2 I 5 



A 



2.s 



™xl/l^i%l- 



(A8) 
(A9) 

(A10) 



a few times of the values given in the Eq. (A7). 

Now we turn attention to the di-Z mode. Compared to the di-phonon mode, the difference 
lies in the additional terms in the polarization vector substitution: e M e* — > g^ v — PiP^/ 771 ^- 
Presumably the resulting change is suppressed by the small parameter m|/m^. (confirmed 
by the h — > Z7 result [39]), then we are justified to ignore this effect at the leading order 
and approximately have 



V(h ->■ ZZ) g 



4^,4 



T(h 77) e*Qt 



4- 



3/2 



c 



(All) 



the phase space suppressing factor is about 0.3 for = 260 GeV, largely it cancels the 
enhancement from coupling ratio. Generally the vector coupling is of the form Vfip^^ipZ^ 



19 



with i)f = Tq — 2Qc sin 2 9w, here Tq is the isospin quantum number of C. 
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